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ABSTRACT 

We revisit the GSC5249-01240 light curve obtained during its occultation by Saturn's North polar 
region. In addition to refractive scintillations, the power spectrum of intensity fluctuations shows an 
enhancement of power between refractive and diffractive regimes. We identify this excess power as due 
to high amplitude spikes in the light curve and suggest that these spikes are due to caustics associated 
with ray crossing situations. The flux variation in individual spikes follows the expected caustic behavior, 
including diffraction fringes which we have observed for the first time in a planetary occultation light 
curve. The presence of caustics in scintillation light curves require an inner scale cut off to the power 
spectrum of underlying density fluctuations associated with turbulence. Another possibility is the pres- 
ence of gravity waves in the atmosphere. While occultation light curves previously showed the existence 
of refractive scintillations, a combination of small projected stellar size and a low relative velocity during 
the event have allowed us to identify caustics in this occultation. This has led us to re-examine previous 
data sets, in which we have also found likely examples of caustics. 

Subject headings: occultations — planets: atmospheres — planets: Saturn 



1. INTRODUCTION 

The occultation of background stellar sources by fore- 
ground solar system objects provides one of the well uti- 
lized probes of distant planetary atmospheres (for a re- 
view, see, Elliot & Olkin 1995). The light curve produced 
during such events allow one to extract detailed informa- 
tion related to the thermal structure of the atmosphere 
with spatial resolutions at the distance of the planet of 
order of a few kilometers. In recent years, increasingly 
sensitive instrumentation has allowed occultation observa- 
tions with spatial accuracies not limited by instrumental 
effects, but rather by physical parameters associated with 
the occultation geometry — such as the distance, the rela- 
tive velocity, and the projected size of the occulted star at 
the planet. 

In Cooray et al. (1998), we presented results from 
an occultation by Saturn's North-polar region of the star 
GSC5249-01240 on November 20, 1995, near the time of 
the solar-ring plane crossing of Saturn. This occultation 
was predicted by Bosh & McDonald (1992) to occur with 
the slowest relative velocity, 0.59 km sec^^, of their 203 
predicted stellar occultations by Saturn between 1990 and 
2000. The small projected diameter of the background 
star, ~ 0.1 km, allowed sub-kilometer spatial resolution of 
the Saturn's atmosphere. 

The light curve of GSC5249-01240 obtained during the 
occultation at NASA's Infrared Telescope Facility (IRTF) 
in Mauna Kea, Hawaii showed an excessive number of 
high amplitude spikes in the light curve. In Cooray et al. 
(1998), we tentatively identified these large amplitude vari- 
ations, of over 100% in some occasions, as due to diffrac- 
tive scintillations. This suggestion was primarily based on 
the fact that the projected size of the star at the distance 
of Saturn, ~ 0.1 km, was much smaller than the Fresnel 
scale associated with these observations, ^ 0.7 km. As the 
light curve probes deeper in to the atmosphere, scintilla- 
tion fluctuations observed in such cases are well described 



through a combination of weak and strong scintillation 
theories based on refraction alone (e.g., Narayan & Hub- 
bard 1988). In the present occultation, due to the smaller 
size of the star when compared to the Fresnel scale, we 
entered a regime for the first time where potentially new 
and interesting signatures may be observed. 

We recommend the reader to Cooray et al. (1998) for 
observational details. There, we also present an analy- 
sis based on the temperature structure of Saturn's North 
polar region. Here, we study the light curve under the con- 
text of scintillation theories based on both diffraction and 
refraction and discuss the power spectrum of intensity fluc- 
tuations as well as individual high amplitude spikes. We 
show that the light curve contains scintillation caustics in- 
volving multiply imaged events; the presence of caustics 
in our light curve agrees with theoretical predictions made 
previously by Goodman et al. (1987) with regards to caus- 
tics in pulsar light curves due to interstellar scintillation. 
Similarly, while not recognized as such, the same caus- 
tics events have also been predicted by French & Lovelace 
(1983) due to ray crossing situations associated with den- 
sity fluctuations generated by gravity waves. Our light 
curve probably provides the best evidence for caustics — 
whatever their cause. 

2. SCINTILLATION CAUSTICS 

The application of scintillations to stellar occultation 
light curves, either due to turbulence or coherent gravity 
waves, has been widely discussed in the literature (e.g.. 
Young 1976; French & Lovelace 1983; Hubbard & Narayan 
1998). In general, scintillation processes are expected from 
fluctuations in the underlying density field, or more appro- 
priately the refractivity in the case of stellar occultation 
light curves, with a power-law spectrum of the form: 



Q{k) - Qo [fc' + fcL] 



-P/2 



exp 



'2fc2 



1 



2 



^OUt k <i kr 



(1) 



where k is the spatial WcLve-number and. /cout 

and kc are 

outer and inner cut-off scales of the fluctuation power spec- 
trum. In the case of an exponential atmosphere, the outer 
cut-ofi^ scale is related to the atmospheric scale height, H 
(French & Lovelace 1983). 

Between the two cut-off scales, the usual assumption of 
a power law in density fluctuations leads to two separate 
and well known regimes of scintillations involving refrac- 
tion and diffraction. If the power density spectrum has an 
inner scale cut off between the refractive and diffractive 
regimes, Goodman et al. (1987) showed that there is ex- 
cess power at the same intermediate regime due to caustic 
events. In the case of stellar occultations, these caustics 
can be associated with ray-crossing situations, either due 
to turbulence or gravity waves (French & Lovelace 1983). 
The flux power spectrum in the presence of caustics can 
be broken to three regimes, refractive, intermediate, and 
diffractive with 
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respectively, where a = /3 — 2. 
Narayan (1985), we write 

r[(4-a)/2] 
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The various length scales involved in Eq. 2 are (1) the 
Fresnel scale, rp = (Ad/27r)^/^, when wavelength of ob- 
servations is A and distance to the scattering screen is 
d, (2) the refractive scale, Trof — dd, when the scatter- 
broadened image size is 9 (Eq. 2.10 of Narayan & Hub- 
bard 1988), (3) the diffractive scale, rdif = rp/vref, that 
determines the fringe separation, (4) the inner cut off 
scale of density fluctuations, Tc = 2Tr /kc, (5) the cut off 
scale of refractive scintillations, rcref = ''"ref if /? < 4 or 
Tcrei = r•,ef(r■c/r,ef)(^-*)/^ ff /? > 4, and, finally, (6) 
the cut off scale associated with diffractive scintillations, 
Tcdif — rp/fc- These length scales are inter-related such 
that riofTdif = ''c^'cdif = Tp. Here, length scales, r, are re- 
lated to wave numbers via fc = 27r/r. In the case of stellar 
occultations, note that the Fresnel scale decreases as one 
probes deep in to the atmosphere. 

In the presence of an inner-scale cutoff, which lies be- 
tween the refractive and diffractive scales, kcrd < kc < 
fccdif, one expects the power spectrum of intensity fluc- 
tuations, P{k), to scale as in the same intermedi- 
ate regime between diffractive and refractive scintillations. 
The presence of refractive scintillations requires that the 
outer cut off scale satisfy the condition rout ^ ^'ref (Good- 
man et al. 1987). 

In Fig. 1, we show the light curve of GSC5249-01240 
during the occultation by Saturn. The existence of high 
amplitude spikes are clearly visible. The relevant parame- 
ters for the GSC5249-01240 event are: relative velocity. 
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Fig. 1. -- The GSC5249-01240 light curve during the occulta- 
tion by Saturn's North polar region and discussed in Cooray et al. 
(1998). The existence of high amplitude short duration spikes are 
clearly visible. Note that these appear as pairs, which is an indica- 
tion for caustics. The labeled spikes are shown in Fig. 3. 



v± = 0.586 km sec~^, distance to Saturn, d = 9.194 
AU, wavelength, A = 2.28 /xm, and integration time, 
tint = 0.947 sec. The Fresnel scale is ^ 0.71 km while 
the resolution of the light curve is ^ 0.55 km. Based on 
the spectrum of the star, the projected radius of the star 
at the distance of Saturn was determined to be of order 
0.1 km. For the rest of the discussion, we renormalize the 
light curve based on model fits presented in Cooray et al. 
(1998) and remove the variation in mean flux and only 
consider fluctuations with respect to this mean. 

The power spectrum of intensity fluctuations, P{k), is 
presented in Fig. 2 where we plot k'^P(k), which represents 
equal weighing in logarithmic wavenumber intervals. In 
general, this power spectrum is expected to show two dis- 
tinct peaks associated with refractive and diffractive scin- 
tillations (see, for e.g.. Fig. 3 of Goodman et al. 1987). 
The data clearly show the existence of the first peak asso- 
ciated with refractive scintillations. In the presence of an 
inner-scale cut off to the underlying density field fluctua- 
tions, we expect excess power beyond the refractive peak 
and a flattening of power at large wave numbers. We see 
an evidence for such an increase in power and flattening as 
shown in the bottom plot of Fig. 2. Here, we have focused 
on only the large wave number regime that is of interest. 

The two solid and dashed curves show fits to the power 
spectrum under the assumption of /3 = 11/3, a Kol- 
mogorov type spectrum, and P — 4 in Eq. 2. With the 
slope fixed, the model fitting involves two unknown free pa- 
rameters related to Qo and kref- In order not to bias model 
fitting with respect to the excess, we make use of the power 
spectrum only out to wave numbers less than 3 km~^; 
While we can obtain a fit to the whole power spectrum, we, 
however, cannot explain the apparent fiattening of power 
beyond the turnover at wave numbers greater than ~ 4.5 
km"-'^. Note that the flattening cannot be ascribed to sim- 
ply noise in the power spectrum; white noise fiuctuations 
will lead to scale independent power with P{k) ~ constant 
such that in a k^P{k) plot, one would see a sharp rise. 



3 



; p^4 

r P=ll/3 




k{km ) 



Fig. 2. — The power spectrum of intensity fluctuations. The 
soUd and dashed curves are based on fits to the power spectrum 
with wavenumbers less than 3 km~^ and using Eq. 2 with assumed 
power-laws for the underlying density-field fiuctuation power spec- 
trum. The peak, at wavenumbers near 1 km~^, is associated with 
refractive scintillations, while an additional peak, which we have not 
probed, is expected at high wavenumbers due to diffractive scintilla- 
tions. In the bottom panel, we concentrate on the high wave-number 
region where, to the right, we see an indication for the flattening of 
power at wavenumbers above 4 km~^. Such a flattening is expected 
when caustics are present in the light curve in addition to simple 
refraction (Goodman et al. 1987). 

Using fits to the power spectrum based on theoretical 
expectations, foUowing Eq. 2, we determined that for a 
P = 11/3 power law, rj-of — ?'crcf — 5.7 ± 1.2 km, and 
Qq — 380 ± 130, at the one-sigma confidence level. For 
a (3 — A power law, these numbers are 5.2 ± 1.1 km and 
370 ± 140 respectively. The errors in fitted parameters are 
estimated by allowing the normalized light curve to vary 
randomly up to ± 5%; this accounts, rather conservatively, 
relative photometric errors and any potential uncertainties 
associated with the light curve. Our discussion is not sub- 
ject to absolute calibration uncertainties, as we are only 
interested in relative changes. The associated diffraction 
scale is rdif ~ 0.09 km or fcdif ~ 70 km~^; we do not probe 
such small length scales, suggesting that we have not ob- 
served regular diffractive scintillations in the present data. 

We can determine the inner cut-off scale by comparing 
the amplitude of flattening to that of the refractive peak 
in the power spectrum. For the /3 — 11/3 and 4 power- 
laws, we determine this cut off to be at rc ~ 4.0 ± 1.5 km 



and 4.4 ±1.6 km, respectively. The caustic regime extends 
from the cut off scale associated with refraction, rcrof ~ 6 
km, to the cut off scale associated with diffraction, Tcdif ~ 
0.12 km. Following Goodman et al. (1987), the expected 
flux variations associated with caustic events are given by 
(AF/F)'^ ~ 2(rc/rrof)^ ln(rcr'i.of/r'|,). Based on derived 
values, we find (AF/F) ~ 1.9 which is consistent, though 
some what higher than, the observed amplitudes of spikes 
in the light curve. 

The presence of caustics can also be verified by inspect- 
ing the intensity profile of individual events. In general, 
individual caustics are expected to show the 1/Vd behav- 
ior in the flux variation with distance, d, as the source 
moves away from the critical line. In the case of stellar 
occultation light curves, such a distance dependence is ex- 
pected in the observer plane such that the intensity varies 
as I{d) = Iq + Q{d)aod~'^/^ , where Q{d) is the step func- 
tion, flo is the caustic strength, and /q ~ is the back- 
ground flux. This intensity dependence, at Saturn, can 
be written using the fact that the mapping between the 
source and observer plane involves d = rcref (a^/^'c)^ and 

do ~ ^c/''cref (Goodman et al. 1987), where x is measured 
distance in the scatterer plane. This leads to an intensity 
decrease for high amplitude spikes as I{x) ^ Aq/x, where 
^0 ^ ^c/''crcf ~ addition to this one-sided drop 

in flux, caustics are also expected to show the presence 
of diffraction with an oscillatory pattern superimposed on 
the intensity variation. 

In Fig. 3, we summarize our model descriptions of indi- 
vidual high amplitude spikes. While the expected width of 
the caustic is a priori known, ^ (rrcffp/riy/^ (Goodman 
et al. 1987), we have allowed for smoothing with a Gaus- 
sian proflle, exp(— a;^/2(T^), and treated the width as an 
unknown parameter. Individual caustic curve model fits 
to spikes are determined by three parameters: the caustic 
strength Aq, smoothing width, and the oscillation spacing 
due to diffraction. From fits to six spikes, we determine 
Aq — 1.8±0.4, consistent with the expectation value based 
on parameters derived from the power spectrum analy- 
sis. The smoothing width, cr, is determined to be of order 
0.5 ± 0.2 km and is consistent with the expected width 
of order 0.45 km. Note that the caustic width is higher 
than the size of the star, ~ 0.1 km, based on its spectrum; 
While Goodman et al. (1997) considered the case of a 
background point source, in our case with the finite size 
for the star, the final caustic width is likely to be smoothed 
by the projected star size at the distance of Saturn. 

Note that caustic spikes appear in pairs in the light 
curve; such a pairing may lead one to conclude the pos- 
sibility that the occultation, in fact, involves a binary in- 
stead of a single star, as was observed in the occultation of 
/3-Scorpii by Jupiter (Elliot et al. 1976). There is one ob- 
servation that rules out this possibility; in our spikes, we 
see the alternating right- and left-handed behavior associ- 
ated with the one-sided decrease in fiux. This is naturally 
associated with caustics. 

Before discussing our main results, we note that previ- 
ous studies on refractive scintillations in occultation light 
curves have shown evidence for anisotropy in underlying 
fluctuations such that in Eq. 1, fc^ = p^k'^+ky with p > 1 
(Narayan & Hubbard 1988). While for the present discus- 
sion we ignored the presence of an anisotropic medium, we 
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plan to address this issue as part of a detailed analysis of 
several light curves of the same occultation event obtained 
with the Hubble Space Telescope. 
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Fig. 3. — Individual spikes as labeled in Fig. 1. The model 
fits, shown as solid lines, are based on caustic descriptions with an 
oscillation due to diffraction. The pairs, A&B and C&D, show the 
expected right and left-handed one-sided decrease in flux associated 
with caustics. 

3. DISCUSSION 

Among the two processes for strong scintillation, refrac- 
tive scintillations have been considered in the past as a 
source of fluctuations in the stellar occultation light curves 
associated with giant planet atmospheres. The second, 
diffractive scintillations, was not expected to contribute to 
intensity fluctuations given the required high resolution to 
see its effects and the small stellar size required not to aver- 
age out the diffractive pattern (e.g., Narayan & Hubbard 
1988). In the case of the occultation of GSC5249-01240 
by Saturn, the projected size of the star was smaller than 
the Fresnel scale and the relative velocity during the event 
was of order 0.5 km sec~^. This, in combination, allowed 
sub-km spatial resolution. 

The light curve contains distinct signatures of sharp and 
high amplitude spikes; such a high number of them has not 
been observed previously. While we explored the regime 
below refractive scintillations, the presence of diffractive 
scintillations require spatial resolutions less than 0.1 km, 
below what we have probed. The diffractive scintillations 
are expected to produce flux variations of order 100% at 
regular intervals. In our case, we see high amplitude spikes 
at irregular intervals with normalized fluxes well over unity 
suggesting that the observed phenomenon is not typical. 

The spikes, however, can be described through caus- 
tics whose presence require an inner cut off to fluctua- 



tions in the underlying density power spectrum (Good- 
man et al. 1987). Such a cut off may also be responsi- 
ble for ray crossing one needs for the formation of caus- 
tics (French & Lovelace 1983). Following the notation of 
French & Lovelace (1983), the present light curve has an 
intensity fluctuation strength, A, of 14 and a normalized 
star size of — O.OSrp-. The strength of intensity fluctu- 
ations suggests fractional fluctuations in the refractivity, 
/S.v/v, of order 4 x 10^^. If underlying fluctuations are 
due to gravitational waves, under geometrical arguments, 
ray crossing is expected when vertical wavelengths are less 
than 27ri?/(Ai^/P)-2/3 - H/Q or 7 km since 7J ~ 40 km 
(Cooray et al. 1998) for this region of the atmosphere. 
This value is the constraint when only a single wave is 
present in the atmosphere. The case with a spectrum of 
waves is analogous to that described by Goodman et al. 
(1987) involving turbulence. 

The presence of caustics is aided by several observa- 
tions that involve the detection of excess power in the 
power spectrum of intensity fluctuations and the expected 
flattening of the same power spectrum in the intermedi- 
ate regime between refractive and diffractive scintillations. 
The intensity variation in individual spikes is also well de- 
scribed through the caustic description with a modulation 
due to diffraction. The parameters from the power spec- 
trum are consistent with model descriptions to the spikes 
aiding the suggestion for the presence of caustics. While 
we have not investigated in detail, we may have found oc- 
casional examples of caustic spikes in previous occultation 
light curves, such as the "camel" spike of Elliot & Veverka 
(1976). These, however, do not contain the signature of 
diffraction fringes as observed in the present occultation. 

From our data, we determined the inner scale cut off of 
the underlying density fluctuation power spectrum to be 
around 4.5 km. The density fluctuations are generally con- 
sidered to be caused by processes involving either turbu- 
lence or gravity waves (e.g. Elliot & Veverka 1976; Young 
1976; French & Lovelace 1983; Narayan & Hubbard 1988). 
It is likely that the cut off may be associated with viscous 
dissipation of wave energy. In Cooray et al. (1998), we 
found an increase in temperature of ~ 15 K between ~ 10 
and 1 /i bar of the region probed by this occultation. Us- 
ing our cut off scale and parameters appropriate for this 
occultation, we extend the calculation by Roques et al. 
(1994) on heating by wave decay to estimate an energy 
transfer rate of ^ 10^^ ergs cm~^ s~^. 

The detection of caustics and, for the first time, diffrac- 
tion fringes in an occultation light curve suggest that we 
entered a regime which was previously unexplored and un- 
expected to be present in occultation hght curves. We be- 
lieve that it was a combination of the small velocity, high 
sampling rate, and the small projected size of the occulted 
star that allowed us to probe the intermediate regime be- 
tween refractive and diffractive scintillations for the first 
time. We suggest events with such favorable conditions be 
given additional priority for observation. 
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